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Abstract. In this paper we review recent discovery of millisecond pulsars (MSPs) in eccentric binary systems. Timing these 
MSPs we were able to estimate (and in one case precisely measure) their masses. These results suggest that, as a class, MSPs 
have a much wider range of masses (1.3 to > 2Mq) than the normal and mildly recycled pulsars found in double neutron 
star (DNS) systems (1.25 < M p < 1.44M Q ). This is very likely to be due to the prolonged accretion episode that is thought 
to be required to form a MSP. The likely existence of massive MSPs makes them a powerful probe for understanding the 
behavior of matter at densities larger than that of the atomic nucleus; in particular, the precise measurement of the mass of 
PSR J1903+0327 (1.67±0.01M©) excludes several "soft" equations of state for dense matter. 
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INTRODUCTION 

In recent years, more than a dozen MSPs have been 
found in binary systems with eccentric orbits. These new 
and unexpected discoveries will allow precise measure- 
ments of the masses of several MSPs. The aim of this 
review is to convey the reasons why this is an exciting 
development. 

In the first two sections, we will provide the context 
for this work. We describe how binary pulsars form; 
this will allow us to understand why they are found in 
two main distinct groups with different spin and orbital 
characteristics. We then describe briefly how we measure 
neutron star (NS) masses (and, in a few cases test general 
relativity, GR) with radio timing. 

In the third section, we enter the core of the review, 
i.e., we describe the recent discovery of many MSPs in 
eccentric binaries. Most of these are located in globular 
clusters. We also review the properties of an intriguing 
new system, PSR J 1903+0327, the first eccentric binary 
MSP in the galactic disk. We finally review what can be 
learned from these new binary systems. 

We assume that the reader is familiar with the concept 
of NSs, and has an idea of what pulsar is: a NS with 
strong, anisotropic electromagnetic emission. For distant 
observers, this emission is modulated by the rotation of 
the object in a repeatable, clock-like fashion, much like 
the light "pulses" of a lighthouse. In what follows, we 
concern ourselves with radio pulsars only. 



FORMATION OF BINARY PULSARS 

Of a total of 1826 rotational -powered radio pulsars listed 
in the ATNF catalog JH, 141 are found in binary systems, 



which makes them relatively rare. Their formation is 
described in Lorimer (2008) and references therein. Fig. 
1 is used there to summarize the two main channels for 
the formation of these systems. 

At the top left, we start with a binary system with 
two main-sequence (MS) stars where at least one of the 
components has > 8M . Such stars are bound to 
end their lives as supernovae (SN) after only a few Myr 
of evolution (next below). This can disrupt the binary 
system, although the exact probabilities depend on the 
orbital parameters, the previous masses of the binary 
components and the magnitude and direction of the SN 
kick. 

If the system survives the first SN, we observe a NS 
orbiting a MS star. There are at least four known ex- 
amples of binaries at this particular evolutionary stage 
where the NS is observed as a radio pulsar. As the com- 
panion evolves, it eventually fills its Roche lobe^ At this 
stage, transfer of matter from the companion to the NS 
creates an "X-ray binary'0. The main consequence for 
the neutron star is spin-up: the accretion of matter from 
the companion transfers angular momentum from the or- 
bit to the NS. This can make the NS reappear as a radio 
pulsar, in which case it is described as a "recycled" pul- 
sar. 



1 The Roche lobe is the region around a member of a binary system 
where matter can still be gravitationally bound to it. 

2 X-ray emission is generated by the hot gas in the accretion disk 
around the NS, in the direct impact of matter with the surface of 
the NS and in some cases in thermonuclear deflagration of matter 
accumulated at its surface, giving X-ray binaries a rich observational 
phenomenology. High-energy emission (in X and gamma-rays) is also 
generated by the interaction of the pulsar wind with the companion's 
wind. 
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FIGURE 1. Formation of millisecond pulsars. From Lorimer (2008). 



The nature of these X-ray binaries and what happens 
during subsequent evolution is determined mainly by the 
mass of the companion. 

1. If the companion's initial mass is also above the 
8M0 threshold, then it is also fated to explode as a 
SN and form a second NS. If the SN kick is cor- 
rectly aligned, the two NSs remain bound, form- 
ing a DNS system. One of the NSs might then 
be detected as a recycled radio pulsar, as in the 
case of PSR B 19 13+16 Jl, and the other might 
be detected as a young radio pulsar, as in the 
case of PSR J1906+0746 [4]. In one known case, 
PSR J0737-3039, both NSs are detectable as radio 
pulsars ]5[]. 

These systems always have eccentric orbits: the 
NSs behave like point masses, so no tidal circular- 
ization happens after the second SN. As we will see 
below, these two features make these systems espe- 
cially useful in the study of gravitation. 

2. If the companion's mass initial mass is below the 
SMq threshold the companion will evolve much 



more slowly and eventually form a white dwarf 
(WD) star ]6j]. Because the companion is an ex- 
tended object for up to several Gyr after the first SN, 
the orbit is very likely to be tidally circularized (un- 
less the separation between components is large). 
No SN or other sudden events occur later on that 
can change this state of affairs. For that reason, the 
vast majority of pulsar - WD systems found in the 
disk of the Galaxy (and all those with orbital peri- 
ods smaller than a few hundreds days) have nearly 
circular orbits. 

The different accretion histories have other important 
implications, which matter directly to the work described 
here: 

Spin Periods. If the companion is light and forms 
a WD its evolution timescale is of the order of several 
Gyr. The long accretion episodes made possible by this 
slow evolution result in NS spin frequencies of several 
hundred Hz, i.e, the NS becomes a "millisecond pulsar" 
(MSP). The first example of this class, PSR B1937+21, 



was discovered in 1982 0]. With a spin frequency of 642 
Hz, this was until 2006 the fastest spinning pulsar known. 

If the companion is massive and forms a NS its evolu- 
tionary timescale is of the order of a few Myr. Given the 
much shorter accretion episode, the resulting spin fre- 
quencies for the accreting NS are one order of magni- 
tude smaller than for MSPs. For that reason, we desig- 
nate them here as "centi-second pulsars" (CSPs). 

Magnetic fields. One apparent consequence of ac- 
cretion is a greatly diminished surface magnetic dipole. 
CSPs, and to a greater extent MSPs have magnetic 
dipoles 3 to 5 orders of magnitude smaller than those 
of normal pulsars. This results in a much smaller brak- 
ing torque and a much longer lifespan (of the order of a 
Hubble time) as radio pulsars. 

Another very important effect of the much smaller 
torque is a much cleaner rotation: effects like "timing 
noise" and "glitches", which make the rotational phase 
of most young pulsars unpredictable, is greatly dimin- 
ished or even absent in CSPs and MSPs. This is fortu- 
nate, because the short spin periods allow very precise 
radio monitoring of the spin phase of these pulsars. This 
monitoring shows that some of the latter objects appear 
to be more stable than atomic clocks; this makes them 
extremely useful astrophysical tools (see below). 

Masses. Another consequence of the different accre- 
tion histories has to do with mass transfer. Since the ac- 
cretion episode necessary for the formation of a MSP is 
much longer than for CSPs, one can reasonably expect 
that the amounts of matter accreted by the former are 
much larger than by the latter. However, it is not clear 
whether this should introduce a systematic difference in 
mass between CSPs and MSPs. The recent results on 
MSP mass measurements (see below) shed some light 
on this matter. 



HOW TO MEASURE THE MASS OF A 
NEUTRON STAR 

As we described above, among radio pulsars recycled 
pulsars have the highest timing precision and the most 
stable rotation. Furthermore, a majority of them are 
found in binary systems, where such features can be 
more profitably employed. 

After discovery of a binary pulsar, we can use Doppler 
variations in the spin period to measure the orbital veloc- 
ity changes along the line-of-sight. The situation is anal- 
ogous to that of spectroscopic binaries, where the chang- 
ing Doppler shift continuously changes the wavelengths 
of the star's spectral lines. Because we can only measure 
changes of velocity along the line of sight our knowledge 
of the binary parameters is incomplete. 



Keplerian Orbits 

In the case of a Keplerian orbit, the measured varia- 
tions of the line-of-sight velocity can lead (after a least- 
squares fit to the observations) to a determination of five 
"Keplerian" parameters: the orbital period (P/,), the semi- 
major axis of the orbit projected along the line of sight (x\ 
for pulsars, this quantity is normally expressed in light 
seconds), the orbital eccentricity (e), the time of pas- 
sage through periastron (7o) and the longitude of peri- 
astron (co). Three unknowns remain undetermined, the 
two component masses (mi, mj) and the orbital inclina- 
tion (z). For most spectroscopic binaries (and most binary 
pulsars) these quantities are simply not available. 

However, there is an equation that links these three 
unknowns, the mass function: 

j,_47Z 2 x 3 _ (ra 2 sin/) 3 ^ 
G Pi (mi +m2) 2 ' 

where G is Newton's gravitational constant. This allows 
an estimate of one parameter (say, m%) from assumptions 
for the values of the other parameters ( mi and i). 

Two more equations are needed to solve these 3 un- 
knowns. No other equation can be obtained from the Ke- 
plerian parameters, and generally no more equations are 
readily available. 



Binary Pulsars 

The distinguishing feature of binary pulsars is that 
we can determine the range directly from the measured 
time of arrival (T.O.A.s) of the radio pulses. Such a 
measurement is impossible with spectroscopic binaries. 
This happens because the radio pulses (unlike spectral 
lines) repeat predictably at a single frequency in the star's 
own reference frame (the neutron star's spin frequency). 

This is extremely important because it allows for an 
astounding gain in precision. As an example, when the 
orbit of PSR J2016+1947 was published d], the orbit 
was estimated from variations in the Doppler shift of the 
spin frequency of the pulsar. The projected size of the 
orbit and the eccentricity determined from this method 
werex= (150.70±0.07) lt-s and e = 0.00128 ±0.00016. 
After the correct rotation count was determined for all 
observations of this pulsar, the range of the pulsar relative 
to the center of mass of the binary could be measured di- 
rectly, with a precision of a few km in each instance. As a 
result, we now obtain jc = (150.7730407 ±0.0000009) lt- 
s and e = 0.001479863 ±0.000000016. This represents a 
gain in precision of about 10 s and 10 4 respectively. This 
is the fundamental reason why binary pulsars are supe- 
rior astrophysical tools. 



For most binary pulsars we are still unable to provide 
any extra equations to help solve the nt\, m.2 and i system, 
despite the unmatched precision provided by the accurate 
ranging. 



Post-Keplerian Effects 

In a few cases, the precision provided by pulsar timing 
and the peculiarities of the system are such that small 
"post-Newtonian" deviations from a Keplerian orbit due 
to the effects of GR become detectable in the times of 
arrival of the pulses. These can be parametrized by five 
quantities, known as "Post-Keplerian (PK) parameters": 
the rate of advance of periastron (ft)), the "Einstein de- 
lay" 7 (due to the larger than average gravitational red- 
shift and special relativistic time dilation near perias- 
tron), the rate of orbital decay due to gravitational radia- 
tion Pf, and two parameters that characterize the effect of 
the gravitational field of the companion on the propaga- 
tion of the pulsar's radio signal (a.k.a. "Shapiro" delay): 
the "range" (r) and "shape" (s). If GR is the correct de- 
scription of gravity, the PK parameters are given by: 
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where m p and m c are the pulsar and companion masses, 



M is the total mass, Tq 
and 



GM /c 3 = 4.925490947^is 



(7) 



The main thing to extract from these equations is that 
they depend on precisely measurable Keplerian param- 
eters, but also on combinations of m p , m c and i. There- 
fore, the measurement of two PK parameters provides, 
together with eq.Q] enough equations to solve for m p , m c 
and i. 

The measurement of additional PK parameters allows 
a test of the self-consistency of general relativity, or 
at least a verification that there are no classical con- 
tributions due to one of the components having a fi- 
nite size. The most famous example of a GR test is 
that carried out with the first binary pulsar ever found, 
PSR B1913+16. The system has an eccentric (e = 0.617) 
and compact (Pi, = 7 h 45 m ) orbit, which means that we 



can measure precisely the longitude of periastron (ft) = 
226.57518(4) t 0on 1986 January 14) and how it changes 
with time (ft) = 4.226607 (7)°yr _1 ). The eccentricity 
causes a non-zero 7, which becomes measurable because 
of the fast precession: 7 = 0.004294(1 )s. Assuming GR, 
i.e., assuming that equations|2]and[3]apply, we obtain for 
the CSP m p = 1 .4408 (3 )M @ and for the younger unrecy- 
cled (and undetectable) NS m c = 1.3873(3)M Q H]. 

These masses allow us to predict the other three PK 
effects for PSR B1913+16 using eqs.|4]-|6l again assum- 
ing that GR applies. The compactness of this system and 
the long timing baseline allowed the measurement of one 
of these PK effects, the orbital decay due to the emission 
of gravitational waves: P h = -2.4211(14) x 10 _12 ss _1 
H. This in perfect agreement with the prediction of eq. 
[4] Apart from testing the self-consistency of GR, this 
measurement demonstrated the existence of gravitational 
waves in the Universe. For this measurement the discov- 
erers, Russel Hulse and Joseph Taylor earned the Nobel 
Prize in Physics in 1993. 

A useful way of visualizing these consistency tests is 
by drawing a mass-mass diagram (Fig. [4] represents the 
constraints derived for PSR J1903+0327). In such a dia- 
gram, m p and m c are the two orthogonal axes. The mea- 
surement of each PK parameter limits the possibilities 
for m p and m c to a small band of the mass-mass diagram. 
In the case of DNS systems, where the components are 
point masses, all of these bands must meet at a single 
point, otherwise GR fails the test. That has not happened 
yet: Since the discovery of PSR B 19 13+16 similar tests 
have been made for other systems (for a review on this 
vast subject see [10]) and GR has passed all of them. 
A system discovered more recently, the "double pulsar" 
(PSR J0737-3039)j|l, allows a total of four tests or GR 
from timing alone lUlll . 



Neutron Star Masses 

Both PSR B1913+16 and PSR J0737-3039A are 
"centi-second pulsars". All precise neutron star mass 
measurements made to date come from CSPs located 
in DNS systems. They range between 1.2489(7)M Q for 
PSR J0737-3039B, a "normal" pulsar that was not recy- 
cled at all JH and 1.44O8(3)M for PSR B 19 13+16 M- 
These values are very close to the Chandrasekhar mass 
(~ 1 .4M ), the upper mass limit for white dwarfs, above 
which these objects become gravitationally unstable. 

Is such a narrow distribution of masses also observed 
for MSPs? From our brief discussion on how they form, 



3 The digit in parenthesis represents the uncertainty in the previous 
digit. 



MSPs could in principle be more massive, owing to the 
much longer accretion episode that formed them, but 
until recently no MSPs had their masses determined to 
less than 10% precision and no precise GR tests have 
been carried out in MSPs-WD systems. 

This state of affairs is somewhat surprising, particu- 
larly considering that the timing precision for MSPs can 
be one or two orders of magnitude better than for CSPs. 
However, as we remarked in the section on the formation 
of binary pulsars, systems containing MSPs form with 
very small orbital eccentricities (2 x 10 -3 to 10~ 7 ). As 
remarked above for PSR B 19 13+ 16, it is the eccentricity 
of the orbit that allows a measurement of cb and 7, and it 
is the combination of these two parameters that generally 
provides the precise mass measurements for the NSs in 
DNS systems. For the circular orbits of MSP- WD bina- 
ries, no such measurements are possible. 

So what about the other PK parameters? They don't re- 
quire an eccentric orbit to be measured, but P/, is greatly 
amplified by a) a compact orbit, b) a large orbital ec- 
centricity and c) a massive companion (eq. |4j. The latter 
requirements mean that even for MSP- WD systems in 
compact orbits (where WD masses are of the order of 
0.1-0.2 Mq and the eccentricities are always very small) 
the orbital decay is is much smaller than in DNS systems. 
This makes its measurement quite difficult, despite the 
much better timing precision of MSPs. Nevertheless, or- 
bital period decays due to GW emission have been mea- 
sured and used to estimate MSP masses (e.g., lfl2ll ). 

How about the Shapiro delay? This is a small timing 
effect (typically, of the order of a few fis in the T.O.A.s) 
which, being proportional to m c , is significantly smaller 
in MSP- WD systems than in CSP-NS systems. However, 
the excellent timing precision of MSPs more than makes 
up for the smaller m c . Nevertheless, this effect seldom 
provides a precise estimate for the pulsar mass in CSPs 
or MSPs. 

This situation is somewhat frustrating: MSPs could in 
principle have higher masses than CSPs and normal pul- 
sars, but despite their excellent timing precision, no pre- 
cise mass measurements can be carried under normal cir- 
cumstances. This results in a lack of GR tests made with 
MSP- WD systems, which is a great pity: the difference in 
binding energies between the components makes MSP- 
WD systems potentially very powerful in detectin g vi - 
olations of the strong equivalence principle (SEP) lUOIl . 
the fundamental physical basis of GR. 

This situation has started to change in recent years, as 
described below. 



ECCENTRIC BINARY MILLISECOND 
PULSARS IN GLOBULAR CLUSTERS 

In globular clusters low-mass X-ray binaries are three 
orders of magnitude more common than in the Galaxctic 
disk. Given the unusual stellar density at the cores of 
globular clusters, it is likely the MS stars are captured 
by old NS lurking near the centers of these clusters. 
Evolution of such MS stars leads to accretion into the 
NS, i.e., large numbers of X-ray binaries III3I1 . These 
eventually become MSPs. 

Since 1987, 140 pulsars have been discovered in glob- 
ular clusters (GCsfj 3/4 of the total were discovered in 
the last ten years. As predicted, the vast majority of this 
pulsar population consists of MSPs, a situation that is 
very different than what is observed in the Galaxy. 

As in the Galaxy, the X-ray binary phase produces to 
low-eccentricity orbits, but in GCs gravitational interac- 
tions with passing stars (or even exchange encounters) 
can, in a few cases, produce binary MSPs with highly 
eccentric orbits 1 14]. The first of the very eccentric (e > 
0.3) binary MSPs to be discovered, PSR J0514-4002A, 
a 5-ms pulsar located in NGC 1851 has e = 0.888 d; 
more so than for any recycled pulsar in the Galaxy. 

Such high eccentricities have already allowed the mea- 
surement of cb for all the eccentric binary MSPs with 
well-known orbits (see TableQ]). If the observed cb is en- 
tirely due to the effects of GR (if, e.g., there are no clas- 
sical contributions due to the finite size of the companion 
star), then this allows an estimate of the total mass of the 
binary M (see eq.|2]i. These values are listed in TableQ] 

As we have seen above, two equations (for the mass 
function / and cb) are not enough to determine the three 
unknowns m p , m c and i. In eccentric orbits, we can 
always measure 7 and obtain the third equation needed to 
solve the system. However, to be able to do so, we must 
wait for the orbit to precess: if there is no precession, the 
effect on the T.O.A.s parametrized by 7 can be absorbed 
by a small re-definition of the size of the orbit, x. The 
orbits of the eccentric MSPs are not as compact as those 
of the DNS systems for which we have measurements 
of 7; this means that the relativistic precession is much 
slower. The end result is that it takes many years of 
timing (more than has elapsed since the discovery of the 
eccentric binary MSPs in TableQ} to measure 7. 

Therefore, the m p -m c -i system cannot yet be solved 
for any of the eccentric binary MSPs in globular clusters. 
Despite this limitation, several important things can be 
learned from a measurement of cb alone, as discussed 
below. 



For an updated list, see http://www2.naic.edu/~pfreire/GCpsr.html 



Early Results 

The first measurement of ft) for a MSP was made in 
2003 for PSR J0024-7204H, a 3.21-ms pulsar in the 
globular cluster 47 Tucanae lltjll . The pulsar is in a 
mildly eccentric (e = 0.07) binary system with an or- 
bital period of 2.35 days [17]. The measured periastron 
advance (ft) = (0.066 ± 0.001)°yr _1 > H) implies, un- 
der the assumption of no classical contributions, M = 
(1.61±O.O4)M . 

As mentioned above, this measurement alone is not 
enough to determine m p , m c and i. However, we can 
combine eq. Q] with M (derived from ft)) to determine a 
minimum companion mass from the condition sin / < 1: 

m e >(/M 2 ) 1/3 - (8) 

A maximum pulsar mass can then be derived from m p = 
M — m c ; the idea is illustrated graphically in Fig. [2] In the 
case of PSR J0024-7204H, this implies m p < 1.52M , 
i.e., the mass of this MSP cannot be much higher than 
the masses we find in DNS systems. The important point 
about this result (which was not even briefly mentioned 
in 1 1 K| ) is that it shows that at least some NSs can be spun 
up to MSP periods with relatively small {dm < O.2M ) 
amounts of matter. 

The discovery of PSR J1909-3744 d] in a search 
of intermediate Galactic latitudes provided a similar re- 
sult. This 2.9-ms pulsar is interesting because it has the 
narrowest pulse profile for any known pulsar (W50 = 
43/is); as a result, it is one of the most precisely timed 
MSPs. Furthermore, the orbital inclination is close to 
90°. This fortunate combination of characteristics al- 
lowed the most precise measurement of a Shapiro delay 
to date. The precise values for r and s (and f) provide un- 
ambiguous values for m c , i and m p = (1 .438 ± O.O24)M 
j20ll . Again, this is consistent with the masses measured 
in DNS systems. Similar results have been obtained for 
other MSPs (see Tableffl)- 

At the start of 2005, it was known that at least some 
MSPs have masses similar to the NSs found in DNS sys- 
tems. Spin-up to MSP periods can definitely be achieved 
with small amounts of matter. 



Recent results 

In 2005, IH discovered 21 new MSPs in the glob- 
ular cluster Terzan 5. Nine other MSPs have been 
found since, including the fastest-spinning pulsar known 
PSR J1748-2426ad, henceforth Ter 5 ad JH], making 
for a total of 33 pulsars in this cluster alone. A total of 
seventeen pulsars are members of binary systems; six of 
these have very eccentric (e > 0.3) orbits. This means 
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FIGURE 2. Mass-Mass diagram for PSR B1516+02B, lo- 
cated in the globular cluster M5. The combination of relativistic 
<X> and the mass function provides an unambiguous maximum 
for m p and minimum for m c . The probability distribution func- 
tions on top and left are calculated assuming a probability for 
the orbital inclinations that is constant in cos; space. 50% of 
the probability lies between the ( = 90° and the i = 60° lines. 



that Terzan 5 alone contains half the known population 
of eccentric binary MSPs. 

At the time of publication, only two eccentric binary 
pulsars Ter 5 I and J had phase-coherent timing solu- 
tions (i.e., an unambiguous pulse count for all detected 
pulses) with precise measurements of ft). Assuming the 
effect to be relativistic (a good assumption, because if 
the companion was extended the orbits of these two pul- 
sars would circularize in about 10 5 yr) the total system 
masses are (2.17 ± 0.02 )M and (2.20 ± O.O4)M . The 
upper mass limits derived for these two pulsars are very 
similar, 1.96M . 

In principle, this is consistent with the masses being 
within the range observed for the NSs in DNS systems. 
However, it is likely that the pulsar masses are signif- 
icantly higher than 1.44 M . The small mass functions 
for these systems imply that the companions are likely to 
contribute little to the total binary mass. 

The argument is probabilistic, and it arises in part from 
assuming that there is no preferred orientation for the bi- 
nary orbits. For randomly aligned orbits, it is much more 
likely that the orbital inclination is close to 90° (edge-on 
orbits) then to 0° (face-on orbits). The reason for this is 
that for face-on orbits only one possible orientation exists 
(orbital plane = plane of the sky), while for edge-on or- 
bits there is an infinite number of possible orbital planes 



TABLE 1. MSP binaries . Notes: a) Binary systems are sorted according to the total estimated mass M. b) Methods are: P), - relativistic orbital decay, r, s - 
Shapiro delay "shape" and "range", "Opt" - optically derived mass ratio, plus mass estimate based on spectrum of companion, 6) - precession of periastron. c) 
This pulsar is not technically a MSP, its spin period is longer than those found in most DNS systems. However, given the similarity of its orbital parameters to 
those of Terzan 5 I, we assume that it had a similar formation history, d) Because of its large companion mass and eccentricity, this system is thought to have 
formed in an exchange interaction flUl . 

Name PSR GC P (ms) P b (days) e f/M @ M/M @ (a) M c /M e M p /M @ Method (b) Ref. (c) 

MSP Mass Measurements 



J0751+1807 - 3.47877 0.26314 0.00000 0.0009674 - - 1.26+[^ P h ,s [12J 

J1911-5958A NGC6752 3.26619 0.83711 <0.00001 0.002688 1-58+jj;^ 0.18(2) l-40lg ^ Opt. [2£] 

J1909-3744 - 2.94711 1.53345 0.00000 0.003122 1.67+| 0.2038(22) 1.438(24) r,s [20] 

J0437-4715 - 5.75745 5.74105 0.00002 0.001243 2.01(20) 0.254(14) 1.76(20) r,s [26] 

J1903+0327 - 2.14991 95.1741 0.43668 0.139607 2.88(9) 1.051(15) 1.74(4) a,s [27] 



Binary systems with indeterminate orbital inclinations 

J0024-7204H 47 Tucanae 3.21034 2.35770 0.07056 0.001927 1.61(4) > 0.164 < 1.52 m [18] 

J1824-2452C M28 4.15828 8.07781 0.84704 0.006553 1.616(7) > 0.260 < 1.367 (b [29J 



J1748-2446I Terzan 5 9.57019 1.328 0.428 0.003658 2.17(2) > 0.24 < 1.96 a [21] 

J1748-2446J (c) Terzan 5 80.3379 1.102 0.350 0.013066 2.20(4) > 0.38 < 1.96 a> [21] 

B1516+02B M5 7.94694 6.85845 0.13784 0.000647 2.29(17) >0.13 < 2.52 m [24J 

J0514-4002A (d) NGC 1851 4.99058 18.7852 0.88798 0.145495 2.453(14) > 0.96 < 1.52 a [30J 

J1748-2021B NGC 6440 16.76013 20.5500 0.57016 0.000227 2.91(25) > 0.11 < 3.3 a [25] 



containing the line of sight. For a system with no known 
orbital inclination, the a priori probability of i\ < i < i% 
is given by cos i\ — cos«2. 

Taking this, CO and its uncertainty into account, we can 
calculate a probability distribution function (p.d.f.) for 
the mass of any MSP and for the mass of its companion. 
For Ter 5 I and J, we obtain median values for the masses 
above 1.7M Q in both cases. Combining the p.d.f.s for 
both pulsars, we obtain a 95% probability that at least 
one of the NSs has a mass over 1.68M . 

Two other results suggest the possibility of even more 
massive pulsars. PSR B1516+02B is a 7.9-ms binary 
MSP in a mildly eccentric (e = 0.14) orbit located in 
the globular cluster M5 B23I1 . After 19 years of Arecibo 
timing the cb for this system has been finally measured 
to good precision l24ll . The total mass of the system 
is (2.29 ±0.17)M Q . This and the small Keplerian mass 
function for the companion imply m p = 2.08^q}qMs 
(these are median and 1-rj limits). There is a 0.7% prob- 
ability that the pulsar has a mass similar to the NSs 
found in DNS systems, and there is a 95% probability 
of m p > 1.72M© (see Fig.0. 

As in the case of Terzan 5 I and J, we have to make the 
assumption that the precession of periastron is relativis- 
tic. In the case of M5B, this is a relatively safe assump- 
tion: HST archival data shows no companion at the loca- 
tion of the pulsar. This implies that it is a small-sized star, 
either a WD or a main-sequence star with m c < O.3M 0! . 
In either case, the contribution of the companion's rota- 
tion to the observed cb is likely to be small |2J]- 

An even more massive NS might be located in the 
globular cluster NGC 6440. PSR J1748-2021B is a 
16.7-ms pulsar with an eccentric (e = 0.57) and wide 
(Pi, = 20.55 days) orbit with a low-mass companion. The 
relatively good precision of its timing and the high ec- 
centricity provide a good measurement of cb B25I1 . Again, 
assuming that there are no classical effects on the or- 
bit, the total mass of the binary is (2.92 ± O.2O)M . To- 
gether with the small mass function measured for this 
system, this implies an extraordinary mass for this pul- 
sar, (2.74±O.2O)M , with a mere 1% probability of 
m p < 2.OM . Until the appropriate optical studies are 
carried out, this result is not as secure as in the case of 
M5B. 
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FIGURE 3. Probability distribution functions (pdfs) for the 
eccentric MSPs binaries in GCs. The mass pdfs of the MSPs in 
the less massive binaries (those with M < 2M ) are represented 
by the dashed curves. The distribution of masses is much 
broader than is found for the components of DNS systems. 
Out of a total of seven systems four seem to be significantly 
more massive than the most massive NS in DNS systems, 
PSRB1913+16. 



that the mass distribution for MSPs is much broader that 
that observed for CSPs and normal NSs in DNS systems 
(see Fig. [3}. About half of the eccentric MSPs in globular 
clusters seem to have large (> 1.6M ) masses. 

It appears therefore that, although some MSPs were 
spun up with relatively modest amounts of matter, that 
was not the case for all MSPs. The reasons for this 
somewhat surprising result are not known at present. 

In the Galaxy, the situation appeared to be different 
until 2008, with PSR J075 1 + 1 807 and PSR J 1 909 - 3744 
indicating relatively low masses and other MSPs yielding 
inconclusive results. The situation has now changed for 
Galactic MSPs as well with the mass measurement of 
PSR J0437-4715 and in particular the discovery of 
PSR J 1 903+0327 JU 



PSR J1903+0327 



MSP Mass Distribution 

The timing results for these MSPs do not prove that 
there are neutron stars with masses significantly in excess 
of 1.44M . It is still possible (but highly unlikely) that 
there are classical contributions to the observed cb for 
most of them. They are, however, highly suggestive: 
plotting the p.d.f.s. and taking them at face value, we see 



PSR J 1903+0327 was a the first MSP found in the on- 
going ALFA pulsar survey ll3~ill . It is a 2.15-ms pulsar in 
a 95-day orbit with a ~ 1M main-sequence star com- 
panion. That is very unusual, but even more unusual is 
the orbital eccentricity of the system: e = 0.44. The stan- 
dard evolutionary scenarios briefly described at the start 
of this review cannot explain the formation of such a 
system; this is one of the reasons why this MSP is so 
interesting. As in the case of the isolated MSPs in the 




FIGURE 4. cosi-m c and m p -m c diagrams for PSR J1903+0327. The contour levels enclose 99.7%, 95.44% and 68.3% of all 
probability. The light contour levels are derived from a y} map calculated from r and s only; the heavy contour levels are calculated 
assuming that the observed 0) is relativistic. Top and Right: projected 1-D p.d.f.s for cos;, m p and m c . The pulsar and companion 
p.d.f.s are much narrower when we take the CO into account, but entirely within the regions predicted using r and i - alone. Orange 
lines: regions of the diagrams consistent with the PK parameters and their 1-CT uncertainties estimated by TEMP02. 



Galactic disk the formation of this system is still not 
well understood [27]. As an example, it is possible that 
PSR J 1903+0327 formed in the same way as the iso- 
lated MSPs, i.e., by somehow eliminating its mass donor, 
which was much closer to the pulsar than its present MS 
star companion. The latter object has been there from 
the start, formerly as the outer element of a hierarchical 
triple; it has never interacted significantly with the MSP. 

The other reason why this binary MSP is so interest- 
ing is that its unusual characteristics mean that 3 PK pa- 
rameters can be measured precisely. The eccentricity al- 
lows a precise measurement of to and the large compan- 
ion mass (m c ~ 1M ) allows a measurement of r and s. 
Combining to and s as measured at the end of 2007 l27ll 
obtained m p = (1.74 ± 0.04)M (;) ; with the ever-present 
qualifier that we are assuming 6a to be relativistic. At this 
time only a small fraction of the orbit had been measured 
at high timing precision (2.2 GHz) with Arecibo, and r 
could not be measured precisely. 

To measure to and s more precisely and verify whether 
to is relativistic, we have started a dense timing campaign 
with Arecibo. The idea is to measure r precisely and see 
if it is consistent with the companion mass we derive 
from to and s. 



When the first 2.2-GHz orbit was completed with 
Arecibo, the pulsar mass estimate decreased to (1.67 ± 
0.01 )M©, but it has been stable at that level the last 
18 months. The companion mass derived from the lat- 
est values of 6a and s is (1.028 ±0.004)M©. The latest 
value measured for r is (1.03 ± 0.04 )M . This agree- 
ment apparently confirms the assumption that the 6a is 
relativistic. However, because of the lower precision in 
the measurement of r, it is impossible to exclude small 
contributions to 6a. At the moment, if we calculate the 
p.d.f. for the mass of PSR J1903+0327 based on r and s 
alone (i.e., assuming nothing about 6a), we obtain m p — 
(1.67±0.11)M© and a 98.4% probability that the mass 
is above 1.44M© (see Fig. @J. The precision of r is still 
improving significantly with continued timing. 

The mass value derived for PSR J 1903+0327 from 
to and s is the most precise MSP mass ever measured. 
The likely contribution from other effects to to is now 
being evaluated, but it is clear by now that is is likely 
to be small. After more than 30 years of searches, this 
is the first precisely measured NS mass larger than that 
of PSR B 19 13+ 16 and the first that is incompatible (and 
significantly above) the Chandrasekharmass. This pulsar 
proves that accretion can significantly increase the mass 



of MSPs compared to the NSs in DNS systems, i.e., the 
mass distribution for MSPs is definitely wider than the 
NS mass distribution in DNSs. 



STUDIES OF SUPER-DENSE MATTER 

Why is a precise measurement of MSP masses so impor- 
tant? First, because of accretion, they can be more mas- 
sive than NSs in DNS systems. If they are, then they test 
our models of how matter behaves at the center of NSs. 

Given their small size (R ~ 10 km) and large mass 
(M ~ IAMq), NSs contain some of the densest matter 
in the Universe: the core can have densities of several 
times 10 14 gcra~ 3 . They are therefore unique astrophysi- 
cal laboratories for testing theories of nuclear matter un- 
der high pressures and densities well in excess of that of 
the atomic nucleus. 

Because matter at such densities are not readily avail- 
able on Earth, its microscopic composition and the re- 
lation between macroscopic quantities like pressure and 
density (the equation of state, or EOS) are essentially un- 
known. As Fig.[5]shows, there are large variations in pre- 
dicted radii and maximum masses for different candidate 
EOSs. These reflect basic uncertainties about the micro- 
scopic behavior and composition of matter at and above 
nuclear density. Some EOSs, such as GS1, assume that 
large percentages of matter are in exotic states (hyperons, 
deconfined quark matter), any of which produces a de- 
crease in the nucleonic (protonic and neutronic) degener- 
acy pressure at any given density, since a larger variety of 
particles are present. This means that matter is relatively 
more compressible; this results in smaller maximum stel- 
lar masses. Other EOSs, such as MSO, assume a larger 
fraction of nucleons and have higher pressures for any 
given density. This means that they predict matter to be 
relatively incompressible and result in larger neutron star 
mass limits. 

As we can see in Fig. [5] the mass measured for 
PSR J 1903+0327 is higher than the maximum possible 
mass predicted by some EOS models still being consid- 
ered in the literature. Such models are excluded by the 
PSR J 1903+0327 mass measurement. 

The mass measurement for PSR J 1903+0327, being 
the largest known with certainty, also sets a quantita- 
tive upper limit to its central density, which then cannot 
be exceeded in any neutron star of lower mass. As suc- 
cessively more massive neutron stars are observed, the 
limiting central density gets smaller. Approximately, the 
limit is p max < 36p s (M /M mM ) 2 where p s = 2/7 x 10 14 
g cm~ 3 is the nuclear saturation density. A large enough 
observed mass (2M ) could rule out the appearance of 
exotic phases in neutron stars, such as deconfined quarks 
i33ll : these are indicated by the SQM lines in Fig. [5] 



PROSPECTS 

In Terzan 5, we now have measurements of 6) for a total 
of seven binary systems. Two more such measurements 
have been obtained for two new eccentric binary pulsars 
in M28, PSR J1824-2452C and D JH]. The increased 
number of mass p.d.f.s will assist the statistical studies 
of the MSP mass distribution. 

For the eccentric binary MSP with the shortest orbital 
period, Ter 5 I, 7 is now becoming detectable. This will 
yield a precise MSP mass measurement for a second 
system (after PSR J 1903+0327) where we expect the 
pulsar to be massive. Whether such a high mass will be 
confirmed or not remains to be seen. 

Meanwhile, HST time has already been allocated 
for the study of the environs of these candidate super- 
massive pulsars. The detection (or not) of their compan- 
ions will be important to address the issue of whether 
the observed 6) is purely relativistic or not. In the case 
of M5B, the non-detection of the companion in archival 
images seems to indicate that this is indeed the case, but 
that is not as clear in the case of NGC 6440B. 

For the latter binary, our simulations indicate that it 
will take about 15 years to determine 7 with any useful 
precision and determine the mass of the pulsar unam- 
biguously, but the potential rewards are immense, partic- 
ularly for the study of the EOS: a mass well in excess of 
2M would exclude most of the EOSs now being pro- 
posed. We expect that advances in instrumentation, with 
resultant improvements in timing precision, will lead to 
significant results well before that. 

The prospects for the study of PSR J 1903+0327 are 
very bright, particularly at optical wavelenghts (no pun 
intended here!). Measurements of the spectral line widths 
of the companion might address once and for all the is- 
sue of whether the companion is rotating fast or not and 
whether it can contribute to the observed advance of pe- 
riastron. If not, then we can consider that the combina- 
tion of 3 PK parameters (ft), r, s) provides one test of GR. 
Measurements of the spectral line shifts will also deter- 
mine the mass ratio precisely, providing one extra con- 
straint in the system, i.e., one extra test of GR. 

Finally, new searches are likely to find several more 
eccentric binary MSPs, particularly in GCs. 
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FIGURE 5. Mass-Radius relation for neutron stars. Each black curve represents a family of neutron stars masses and radii 
according to a given equation of state. The region bounded by the Schwarzschild condition R < 2GM/c 2 is excluded by general 
relativity, and that bounded by R < 3GM/c 2 (labeled "causality") is excluded by requiring the speed of sound inside the star to 
be smaller than the speed of light. The mass-shedding limit for the fastest spinning radio pulsar (PSR J1748— 2446ad, with spin 
frequency 716 Hz) is labeled "rotation"; points in the green region below this line are not allowed for that particular pulsar. Stricter 
constraints may arise from X-ray sources like XTE J1739— 285 (dashed curves calculated under different neutron star models) if 
their spin frequencies are confirmed to be higher than that of PSR J748— 2446ad, potentially excluding some equations of state 
(such as GM3) which lie almost entirely below the rotation curve. A recent, precise millisecond pulsar (MSP) mass measurement 
(for PSR J1903+0327) excludes the "softest" EOSs (red horizontal line). Adapted from jH. 
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